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Abstract 

We propose a unified framework for equity and credit risk modeling, where the default time is 
a doubly stochastic random time with intensity driven by an underlying affine factor process. This 
approach allows for flexible interactions between the defaultable stock price, its stochastic volatil- 
ity and the default intensity, while maintaining full analytical tractability. We characterise all risk- 
neutral measures which preserve the affine structure of the model and show that risk management 
as well as pricing problems can be dealt with efficiently by shifting to suitable survival measures. 
As an example, we consider a jump-to-default extension of the Heston stochastic volatility model. 

Key words: default risk, affine processes, stochastic volatility, market price of risk, change of mea- 
sure, jump-to-default 

1 Introduction 

The last few years have witnessed an increasing popularity of hybrid equity /credit risk models, see e.g. 
the recent papers [|2l[51|6l[8l|9l[T0l[T4l[T5l. One of the most appealing features of such models is their 
capability to link the stochastic behavior of the stock price (and of its volatility) with the randomness 
of the default event and, hence, with the level of credit spreads. The relation between equity and credit 
risk is supported by strong empirical evidence (we refer the reader to [61 [T5]| for good overviews of 
the related literature) and several studies document significant relationships between stock volatility 
and credit spreads of corporate bonds and Credit Default Swaps (ll4l[T7l). 

In this paper, we propose a general framework for the joint modeling of equity and credit risk 
which allows for a flexible dependence between stock price, stochastic volatility, default intensity and 
interest rate. The proposed framework is fully analytically tractable, since it relies on the powerful 
technology of affine processes, and nests several stochastic volatility models proposed in the literature, 
thereby extending their scope to a defaultable setting. Furthermore, a distinguishing feature of our 
approach is that, unlike the models proposed in [|2l[6l[El|9l[Hl[l5]|, we jointly consider both physical 
and risk-neutral probability measures, ensuring that the analytical tractability is preserved under a 
change of measure, while at the same time avoiding unnecessarily restrictive specifications of the risk 
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premia. This enables us to efficiently compute several quantities which are of interest in view of both 
risk management and pricing applications. 

The paper is structured as follows. Sect. [2] introduces the modeling framework, while Sect. [3] 
gives a characterisation of the family of risk-neutral measures which preserve the affine structure of 
the model. In Sect.s [4] and [5} we show how most quantities of interest for risk management and 
pricing applications, respectively, can be efficiently computed (we refer the reader to Sect. 2.5 of 
[|29l for more detailed proofs of the results of Sect.s |4jj5]). Sect. [6] illustrates the main features of 
the proposed approach within a simple example, which corresponds to a defaultable extension of the 
Heston lT30ll model. Finally, Sect. [TJconcludes. 



2 The modeling framework 

This section presents the mathematical structure of the modeling framework. Let (Cl, Q, P) be a 
reference probability space, with P denoting the physical/statistical probability measure (we want to 
emphasise that our framework will be entirely formulated with respect to the physical measure P). 
Let T E (0, oo) be a fixed time horizon and W = (W t )o<t<T an Revalued Brownian motion on 

Q, P), with d > 2, and denote by F = (J r t)o<t<T its P-augmented natural filtration. 

We focus our attention on a single defaultable firrrj^J whose default time r : Q — > [0, T] U {oo} 
is supposed to be a (P,¥)-doubly stochastic random time, in the sense of Def. 9.11 of [|32l . This 
means that there exists a strictly positive process A p = (Af) <t<r such that, for all t E [0, T], we 
have P{t > t \ Ft) = P{t > t \ Tt) = exp(— f^X^du). We call the process A p the P -intensity 
of r, thus emphasising the role of the reference measure P. Let the filtration G = (Gt)o<t<T be the 
progressive enlargemenj^of F with respect to r, i.e., Q t := V\ s>t {Fs V cr(r A s)}, for all t E [0, T], 
and let Q = Q T . It is well-known that G is the smallest filtration (satisfying the usual conditions) 
which makes r a G-stopping time and contains F, in the sense that T t C Gt for all t E [0, T). 

The price at time t E [0, T] of one share issued by the defaultable firm is denoted by S t . We 
assume that the G-adapted process S = (S t )o<t<T is continuous and strictly positive on the stochastic 
interval [0, rjand satisfies S'lji-.T] — 0. This means that S drops to zero as soon as the default event 
occurs and remains thereafter frozen at that level. By relying on the Sect. 5.1 of together with 
the fact that all F-martingales are continuous, it can be proved that there exists a continuous strictly 
positive F-adapted process S = (S t )o<t<T such that S t = l{ T >t} St holds for all t E [0, T). We shall 
refer to the process S as the pre-default value of S. 

The pre-default value S is assumed to be influenced by the F-adapted stochastic volatility process 
v — ( v t)o<t<T an d by an IR d ~ 2 -valued F-adapted factor process X = (X f ) <t<T- The process X 
can include macro-economic covariates describing the state of the economy as well as firm-specific 
and latent variables, as considered for instance in ll2~8l . Let us define the process L = (L t )o<t<T by 
L t := log S t and the Revalued F-adapted process V = {V t ) <t<T by V t := (v t ,Xj,L t ) T , with T 
denoting transposition. 

The processes v, X and L are jointly specified through the following square-root-type SDE for the 
process V on the state space IR™ + x ]R d_m , where we let M.™ + := {x E M. m : Xi > 0, Vi = 1, . . . , m}, 

'The present modeling framework can be easily extended to the case of M > 1 defaultable firms if we suppose that 
the random default times {ti, . . . , tm} are F-conditionally independent (see e.g. Il32l . Sect. 9.6). 

2 Due to Lemma 6.1.1 and Lemma 6.1.2 of 0, the fact that P(t > t \ T T ) = P(t > t \ T t ), for all t € [0, T], implies 
that all (P, F)-martingales are also (P, G)-martingales. In particular, W = (Wt)o<t<T is a Brownian motion with respect 
to both F and G. This important fact will be used in the following without further mention. 
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for some fixed m G {1, . . . , d — 1}: 

dV t = (AV t + b)dt + Y,jR t dW t 



V = (v ,Xj,\ogSo) T = v G 



i™ x R d - m 



(2.1) 



odxd x fljd x jg>c(xd anc j ^ j s a (ji a g 0na i ^ x rf) -matrix with elements given by 



and/3 := (&,...,&) G 



ocixc 



where (A, 6, S) G 

= cii + /3jV t , for all t G [0, T], with a := (a l9 . . . , a d ) 
In order to guarantee the existence of a strong solution to the SDE ( |2.1| ), we introduce the following 
Assumption. 

Assumption 2.1. Lef m G {1, . . . , d — 1} and define I := {1, . . . , m}, J := {m + 1, . . . , d} and 

D :— I U J — {1, . . . , d}. The parameters A, b, S, a, (3 satisfy the following conditions: 

(i) h > (Sjj) 2 /3j )i /2/or alii G I; 

( ii) Aij = Ofor all i G I and j G J and A it j > Ofor all i,j G I wzY/z i 7^ j; 
fzzz'J Sjj = Ofor alii G / and j E D with j 7^ z; 

( z'vj = Ofor all i E D and j G J, (3 iti > Ofor all i E I and /5y = 0/or a// i,jEl with i 7^ j; 

(V) «j = Ofor alii G I and aj > — Yl^Li f or a ^ j e ^- 



For any v G 



x 



id—m 



, Assumption 2.1 ensures the existence of a unique strong solution 



V = (V t ) <t<T to the SDE (2.1) on the filtered probability space (fi,£?,F,P) such that V 



and V t G x IR d ~ m P-a.s. for all i G [0, T]. Indeed, the same arguments used in the proof of 
Lemma 10.6 of 11271 give the existence of a unique strong solution V = (V t )o<t<T on IR™ x M d_m , 
while Lemma A. 3 of [|24| together with Ex. 10.12 of [27J implies that V actually takes values in 

this also implies that the matrix R t is 



X 



Yid—m 



. Due to conditions (iv)-(v) of Assumption 



2.1 



positive definite for all t G [0, T\. We shall always assume that Assumption 2.1 is satisfied without 
further mention. 



Remark 2.2. The parameter restrictions imposed by Assumption 2.1 bear resemblance to the canon- 
ical representation of [fT9ll . However, we do not require the matrix £ to be diagonal, since this may 
lead to unnecessary restrictions on the model if 2 < m < d — 2, as pointed out in lfT2ll . 

The following Proposition describes the dynamics of the defaultable stock price process S. 
Proposition 2.3. The process S = (S t )o<t<T satisfies the following SDE on (fi, Q, G, P): 

(d-2 m-1 \ 

s + mi log S t ^ + fi 2 v t + ViXl + Vi x t I dt 
i=i i=i / (22) 

+ 5 t _ a^c^ 1 + 5 t _ ^ Z^^RpdWt - S t _ rfl {r < 4} 

i=2 

wzY/z convention St- log = on {r < £} and where 

1 d 1 1 d 

s:=b d +- (E d>jfc ) 2 o; fc , /xi := /i 2 := v4 d> i + - (S dj i) 2 /?n + - ^ (S^) 2 

fe=m+l 

1 1 ^ 



fc=m+l 



:,»+!> 



(7 



+l,fe 



fc=m+l 
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Proof. Observe first that dS t = l{ T> t-}S t - (dL t + d(L) t /2) — St- dl{ T < t }, due to Ito's formula and 
integration by parts. Equation ( |2.2[ ) then follows from ( |2.1| ) together with Assumption |2. 1 1 by means 



of simple computations. □ 

Remark 2.4. The defaultable price process S has a rich structure, influenced by the factor process 
X in both the drift and diffusion terms. Furthermore, there are three levels of dependence between S 
and the stochastic volatility v: (1) a direct interaction, since v explicitly appears in the dynamics of 
S; (2) a "semi-direct" interaction, since the Brownian motion W 1 driving the process v is also one of 
the drivers of S; (3) an indirect interaction, since 5* and v both depend on the factor process X. 

To complete the description of the modeling framework, we specify as follows the P-intensity 
process A p = (Af )o<t<T and the risk-free interest rate process r = (r t )o< t <T' 

X P :=X P + (A p ) T V t r t :=f + T T V t forallte[0,T] (2.3) 

where A p ,f e R + and A p , T e x {0} d ~ m satisfying the conditions A p + Y.T=i A f > and 
f + Y^iLi > 0- This ensures that the P-intensity and the risk-free rate are correlated and strictly 
positive, since is an unattainable boundary for V\ Vz G I. Furthermore, the linear structure ( |2.3| ) 



permits to obtain analytically tractable formulae for several quantities of interest, as shown in Sect.s 
Wjpj The specification (|2.3[) allows for a direct dependence of A p on the stochastic volatility v, this 



feature being consistent with several empirical observations (see e.g. flUHTl). Observe also that the 
defaultable price process S and the P-intensity A p are linked through the common factor process X. 
Finally, we want to note that the proposed modeling framework generalises to a defaultable setting 
several stochastic volatility models considered in the literature. For instance, defaultable versions 
of the models considered in Q3 and Sect. 4.3 of [|25l can be easily recovered within our general 
setting. 



3 Equivalent changes of measure which preserve the affine struc- 
ture 

The modeling framework introduced in Sect. [2]has been formulated entirely with respect to the physi- 
cal probability measure P. However, since we aim at dealing with pricing as well as risk management 
applications, we need to study the structure of the model under a suitable risk-neutral probability mea- 
sure, formally defined as a probability measure Q ~ P on (O, Q) such that the discounted defaultable 
price process exp(— J Q r u du) S is a (Q, G)-local martingale]^ 

It is important to be aware of the fact that most of the appealing features of the framework de- 
scribed in Sect. [2] may be lost after a change of measure. Aiming at a model which is analytically 
tractable under both the physical and a risk-neutral measure, we shall consider all risk-neutral mea- 
sures Q which preserve the affine structure of(V, r), in the sense of the following Definition. 

Definition 3.1. Let Q be a probability measure on Q) with Q ~ P. We say that Q preserves the 
affine structure of (V, r) if the following hold: 



(i) the process V = {V t ) < t <T satisfies an SDE of the type ( |2.1| ) on (Q, Q, F, Q) with respect to an 



l d -valued (Q, ¥)-Brownian motion W Q = (W®)o< t <T and for some parameters A Q , b Q , S, a, 



(3 satisfying Assumption 2.1 



3 Due to the fundamental result of [20 1, this is equivalent to the validity of No Free Lunch with Vanishing Risk (NFLVR) 
condition for the financial market (S, G), being the process exp(— J r u du) S locally bounded. 
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(ii) the default time r is a (Q, ¥)-doubly stochastic random time with Q-intensity X Q = (\f)o<t<r 
of the form X? = X Q + (A Q ) T V t ,forX Q E R+ and A Q E M™ x {0} d - m with X Q + J2T=i A ? > °- 

We denote by Q the family of all risk-neutral measures which preserve the affine structure of 
(V,r), in the sense of Def. 3.1 The following Theorem gives a complete characterisation of the 
family Q. This result follows from a more general one in Chapt. 2 of 11291 , but we outline a self- 
contained proof for the convenience of the reader. We denote by £ the stochastic exponential and by 

M = (M t 



t 0<t<T 



the (P, G) -martingale defined by M t := l{ T <t}— J T X p du (see e.g. Q, Prop. 5.1.3) 



Theorem 3.2. Let Q be a probability measure on (Q, Q). Then we have Q E Q if and only if 



dQ 
dP 



£ 9dW + jdM 



exp 



. i=l 







T o\dwi- l -j2 

i=l 



fdt 



V 2 



tAT 



JtKdt 




(3.1) 



1{t<t}7t 



where 9 = {9 t )o<t<T and 7 = ipft)o<t<T are F-adapted processes of the following form: 



9 t = 9(v t ) -.= R; 1/2 (9 + ev t ) 



It = l{V t ) 



(X Q — X p ) + (A®-A p ) T V t 
X p + (A P ) T \4 



(3.2) 



for some 9 e R d and 6 G R dxd such that: 

(0 Yl=i ^i,k9k > (Sj,j) 2 /3i,j/2 - bi for alii E I; 

(ii) Ylt=i ^i,k®k,j = 0, for all i E I and j E J, and 2~^fe=i — ~Aij, for all i, j E I with 

for some X Q E R+ and A Q E M.™ x {0} d - m with X Q + YZi A * Q > and if the following equality 



holds P-a.s. on {r > t}, using the notations introduced in Prop. 2.3 

d 



s + fii log St- + ( /i 2 + a 



01 



v t - 



d-2 
i=l 



m—1 



r t + Af(l + 7 t) (3.3) 



i=2 



Proof. Let 9 = (9 t ) < t < T and 7 = ('ft)o<t<T be two F-adapted processes satisfying ( |3.2[ ). Since 9 and 
7 are continuous functions of V and the process V is continuous, hence locally bounded, the process 
Z := £^J 9 dW + J jdMj is well-defined as a strictly positive (P, G)-local martingale and, as a 
consequence of Fatou's lemma, it is also a (P, G)-supermartingale. Moreover, Thm. 2.4 and Remark 
2.5 of lfT3l allow to conclude that E[Zt] = 1, thus implying that Z is a uniformly integrable (P, G)- 
martingale. So, we can define a probability measure Q on (fi, Q) via ( |3.1| ). Part (i) of Def. 3.1 then 
follows from Girsanov's theorem together with ( |3.2[ ), while part (ii) follows from Thm. 6.3 of |[T6ll . 
Girsanov's theorem together with ( |3.2| ) and Prop. 6.2.2 of (3). Finally, the (Q, G)-local martingale 
property of exp(— J Q r u duj S easily follows from Girsanov's theorem together with Prop. 



2.3 



and 

ows 



p.3[ ). Conversely, suppose that Q E Q. The existence of a representation of the form ( |3.1| ) fol 
from Corollary 5.2.4 of 0, while p.2[ ) and ( |3.3[ ) follow from Girsanov's theorem together with Def. 
3.1|and Prop. |2.3[ respectively. □ 



Remark 3.3. The processes 9 = (9 t )o<t<T and 7 = (7t)o<t<T admit the financial interpretation of 
risk premia (or market prices of risk) associated to the randomness generated by the Brownian motion 
W and by the random default time r, respectively. More specifically: 
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(a) The process 9 = (9 t )o<t<T represents the risk premium associated to the diffusive risk generated 
by the Brownian motion W. Since the stock price, its stochastic volatility, the default intensity 
and the interest rate all depend on W through V, the risk premium 9 can be considered as a 
market- wide non-diversifiable risk premium^] 

(b) The process 7 = (7*)o<t<r represents the risk premium associated to the default event or, more 
precisely, the risk premium associated to the idiosyncratic component of the risk generated by 
the occurrence of the default event (to this effect, see also [|51|26l and Sect. 9.3 of 11321 ). 

The importance of explicitly distinguishing between 9 and 7 has been demonstrated in ETTl . Assuming 
7 = means that the idiosyncratic component of default risk can be diversified away in the market, as 
explained in [13111 , and, therefore, market participants do not require a compensation for it. However, 
the jump-type risk premium can be significant when it is difficult to hedge the risk associated with the 



timing of the default event of a given firm. Note that, as can be seen from ( |3.2[ ), the risk premia 9 and 
7 both depend on the common driving process V. 



Due to Thm. 3.2 our modeling framework enjoys full analytical tractability under both the phys- 
ical measure P and any risk-neutral measure Q E Q, thus enabling us to efficiently solve risk man- 
agement as well as a pricing problems, as we are going to show in Sects [4j|5] We close this Section 



with the following fundamental result, which follows from Thm. 10.4 of E71 together with part (i) 



of Definition 3.1 ( |2.1[ ) and Assumption 2.1 For z E C d we denote by ^t(z) and ^s(z) the real and 
imaginary parts of z, respectively, and C™ := {z E C m : 3t(z) E W?\. For Q E Q U {P}, we denote 
by E Q the (conditional) expectation operator under the measure Q. 

Proposition 3.4. For every Q E Q U {P} and for all z E C m x iM. d ~ m , there exists a unique solution 
z), z)) : [0, T\ -> C x C d to the following system ofRiccati ODEs: 



d t <S> Q (t } z) = (bQ) T ^(t,z)+ 1 - p T * Q (t,z)}ta k -\ Q -f^ P 
),z) = 



2 

fe=m+l 



d d 



k=l k=m+l \fi E I 

d 

k=m+l 

ty Q (0,z) = z 

(3.4) 

Furthermore, for any Q E Q U {P}, the following holds for all < t < u < T and for all 

z E C m x iR d - m : 



E Q 



exp(^J\\f + r s lQ^p)ds + z T V,}j T t = exp ($ Q (u - t,z) + * Q (u - t, z) T k) (3.5) 



4 In the context of default-free term structure modeling, in ifTTI the authors demonstrate that the specification ( |3.2[ > has 
a considerably better fit to market data than the simpler market price of risk specifications traditionally considered in the 
literature (see e.g. iflOl [1911221123113011. 
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4 Risk management applications 



Many quantities of interest in view of risk management applications can be computed as conditional 
expectations under the physical measure P. As a first and basic application, let us compute the Q t - 
conditional survival probability of the defaultable firm up to the final horizon T. We denote by $ p (-, •) 
and ^ p (-, ■) the solutions to the Riccati ODEs ([374]) with Q = P. 



Proposition 4.1. For any t G [0,T], the following holds: 

P(r>T\g t ) = l {T>t} exp ($ P (T -t,0) + V P (T - t, 0) T V t ) (4.1) 



Proof. Corollary 5.1.1 of implies that P(r >T\ Q t ) = l {T>t} E[exp(- j^Xfds) | T^. The result 



then follows by applying (3.5 ) with Q = P, z = and u = T. □ 



As can be easily checked from ( |3.4| ), the right-hand side of ( |4.1[ ) only depends on {V 1 : % G /}, i.e., 
on the components of the process V on which the P-intensity A p depends. For computing conditional 
expectations (under the measure P) of more general quantities needed for risk management purposes, 
it turns out to be convenient to introduce the T -survival measure P T ~ P on (Q,G) defined by 
dP T /dP := exp(- J T Af dt) jE [exp(- ^Xfdt)] . 

Lemma 4.2. For any random variable F G L}{P, J 7 ^) and for any t G [0, T] the following holds: 

E [F 1 {T>T} | Qt] =P(r>T\ Q t )E pT [F \ JF t ] (4.2) 

Proof Corollary 5.1.1 of implies that E [Fl {r>T} |^ t ] = l {r>t] E[F exp(- j^Xfds) \T t ]. 
Equation ( |4.2[ ) then follows by using the definition of the measure P T together with the conditional 



version of Bayes' formula (see e.g. [|27l . Ex. 4.9). □ 



Lemma 4.2 shows that the computation of the (^-conditional expectation of an J^-measurable 
random variable F in the case of survival until time T reduces to the computation of the J^-conditional 
expectation of F under the T-survival measure P T , the term P(r > T\Gt) being computed as in ( |4.1[ ). 



Furthermore, the T t -conditional characteristic function of the vector Vr under the T-survival measure 
P T can be computed in closed form, as shown in the following Lemma. 

Lemma 4.3. For any z G iM d and for any t G [0, T] the following holds: 



(pf(z) := E pT [e zTVr \F t ] = exp U p (T-t, z)-<S> p (T-t, 0) + (V p (T-t, z)-V p (T-t, 0))V 4 ) 

(4.3) 



Proof. The definition of the measure P T together with the conditional version of Bayes' formula 
gives E pT [e ^ | F t ] = E [exp (- tf\ P ds + z T V T ) | JFt[ /E [exp (- f?\?ds) | JF t ] . By applying 
(|3.5[) with Q = P, u = T and z G iR d (z = 0, resp.) to the numerator (to the denominator, resp.), we 



then obtain equation (4.3 ). □ 



Due to Lemma 4.2 and Lemma 4.3, we can compute the ^-conditional expectation (under the 
physical probability measure P) of arbitrary functions of the random vector Vr in the case of survival 
by relying on well-known Fourier inversion techniques. As an example, we can explicitly compute 
quantiles of the ^-conditional distribution of the defaultable price Sr in the case of survival. This is 
crucial for the computation of Value-at-Risk and related risk measures. 
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Proposition 4.4. For any x G (0, oo) and for any t G [0, T] the following holds: 



P{St < x,r > T\g t ) = P (r > T \Q t ) 



1 



7T 



oo cWg-iylogz m P 



iff (0,...,0,zy)) 



2/ 



(4.4) 



where P(r > T\Gt) an d yf T (-) are explicitly given in ( |4.1[ ) and Lemma 4.3 respectively. 
Proof. Note that: 

P (<Sr < x, r > T|&) = P (Lt < logx, r > T\Q t ) = P (r > T|&) P T (L T < logx^ 



where the second equality follows from Lemma 4.2 Equation ( |4.4[ ) then follows from standard 
Fourier inversion techniques (see e.g. 11291 . Prop. 2.5.12, and 113311 . Sect. 1.2.6). □ 



5 Valuation of default-sensitive payoffs and defaultable options 

Throughout this section, we fix an element Q G Q. For the purpose of valuing default-sensitive 
payoffs, the u-survival risk-neutral measure Q u , for u G [0, T\, turns out to be quite useful. The 
measure Q u is defined by dQ u /dQ = exp(— j^(r s + \® ) ds) /E Q [exp(— J^(r s + A^) ds)] . For 
u = T, the measure Q T bears resemblance to the T -survival measure P T introduced in Sect. [4j 
except that Q T is defined with respect to some Q G Q and the density dQ T /dQ also involves the 
risk-free interest rate besides the Q-intensity \ Q (compare also with 0, Def. 15.2.2). Following 
the same logic of Sect. |4} we show that many pricing problems can be simplified by shifting to the 
measure Q u , for some u G [0, T). As a preliminary, let us compute the arbitrage-free price LT(t, T) of 
a zero-coupon defaultable bond. We denote by •) and •) the solutions to the Riccati ODE 



system (3.4). 



Lemma 5.1. For any t G [0,T] the following holds: 

T) = l {T>t} exp ($ Q (T - t, 0) + ^ Q (T - t, 0) T V t ) 

Proof. Note first that: 

II(t,T) = exp(— ^ r s c/s)l{ T>r }|^] = l{ r>t }i? Q [exp(— y (r s + Ajp) ds)!^] 



(5.1) 



where the second equality follows from Thm. 9.23 of 11321 . Equation ( |5.1[ ) then follows from Prop. 

□ 



13741 with u = T and z = 0. 

Of course, coupon-bearing corporate bonds can be valued as linear combinations of zero-coupon 
defaultable bonds (see Sect. 1.1.5). More generally, most default- sensitive payoffs can be de- 
composed into linear combinations of zero-recovery and pure recovery payments, the latter being 
paid only in the case of default, see e.g. Sect. 9.4 of [|32l . The next Proposition provides general 
valuation formulas for zero-recovery and pure recovery payments. 



Proposition 5.2. For any t G [0,T] and for any measurable function G 
following hold: 



x 



od—m 



E Q 



-ft r » ds / 



G(V T )l {T>T} \g t =n(t,T)E® [G(V T ) I Tt. 



E Q 



■ ft r ' ds 



G(v T )i {T < T} \g t 



H(u,T)E^[X^G(V u )\^ t ]du 



-> R + the 
(5.2) 
(5.3) 
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Proof. Note first that, due to Thm. 9.23 of 11321 . we can write: 



E Q 
E Q 



■f^r s ds 



G(V T )l {T>T} \g t =l {T>t} E Q 



e-K r * ds G(V T )l {T < T} \g t ] = l {T>t} E Q 



-J t T (r 3 +X?)ds ( 



i '•• • C/( l v ) 7, 



f:(rs+X?)ds X Q G( y u)du 



Equations ( ]5.2[ )-( [5T3] ) then follow by using the definition of the measure Q u , for u G [t, T], together 
with the conditional version of Bayes' formula and also, for ( |5.3| ), with Tonelli's theorem. □ 

For < t < u < T, the following Lemma gives an explicit expression for the .^-conditional 
characteristic function (f{ of the random vector V u under the w-survival risk-neutral measure Q u . Its 
proof follows from ( [33 ) and, being analogous to that of Lemma |4.3[ is omitted. 

Lemma 5.3. For any < t < u < T and for any z G iM. d the following holds: 

ipf(z) := E QU [e zTv - 1 Ft] = exp U Q (u -t,z)- ® Q (u - 1, 0) + (V Q (u -t,z)- V Q (u - 1, 0)) T V t ) 

(5.4) 



By combining Prop. 5.2 with Lemma 5.1 and Lemma 5.3 and using well-known Fourier inversion 
techniques, we can obtain semi-explicit formulas for a wide range of default- sensitive payoffs. In 
particular, we now derive valuation formulas for Call and Put options (issued by a default-free third 
party) written on the defaultable stock S. We denote by H r f(t, T) := E Q [exp(— f r s ds) \ Q t ] = 
E Q [exp(— J^r s ds) \ Q t ] the arbitrage-free price at time t G [0, T] of a zero-coupon default-free bond. 



Corollary 5.4. For any t G [0, T] and for any strike price K > the following hold: 



C K (t,T) := E Q \e-£ r ° ds (Sr- K)+\ Q t 



2tt 



$ {0,...,0,w + iu) 



j£—(w—l+iu) 



(w + iu) (w — 1 + iu) 



du 



P K (t,T) := E Q \e-^ r ° ds (K - Sr) + \ Q t 



K(U rf (t,T)-Il(t,T)) + 



n(t,T) 

2tt 



+oo 



(pf\0,...,0,y + iu) 



j^-(y-l+iu) 



(y + iu)(y — 1 + iu) 



(5.5) 
(5.6) 

du 



(5.7) 



for some w > 1 and y < such that the system of Riccati ODEs ( |3.4| ) has a unique solution for the 
initial conditions z = (0, . . . , 0, w) T and z — (0, . . . , 0, y) T . 



Proof. Observe first that: 



E Q 



e-^ ds (Sr-K) + \Q t 



E Q 



■ St r = ds 



(Sr-K) + l {T>T} \g t =U(t,T)E QT (e LT -K) + \F t 



where the second equality follows from ( |5.2[ ). As in Q and [|27ll . Lemma 10.2, it can be shown that: 

J£— (w— 1+iu) 



(e* - K)~ 



1 

27 



du 



, (w + iu) (w — 1 + iu) 

for some w > 1. Equation ( |5.6| ) then follows by Fubini's theorem (see Cor. 2.5.21 of [|29l for more 
details). Equation ( |5.7[ ) follows by an analogous computation once we observe that: 



E Q 



e-^ ds (K-Sr) + \Q t 



E Q 



e I t T rsds (K-fy + i {r>r} | g t j +KE Q [ e - f?r.d. ^ _ 1{t>t} ) | q% 

□ 
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6 An example: the Heston with jump-to-default model 



In this section, we illustrate some of the essential features of the proposed modeling framework within 
a simple example, which corresponds to a generalisation of the stochastic volatility model introduced 
by Heston [|30l . here extended by allowing the stock price process to be killed by a jump-to-default 
event, in the spirit of [|9]|. 



6.1 The model 

Using the notations introduced in Sect. [2} we let d = 3 and consider the following specification: 



.4 



-k 


-1/2 





(a 















Rt 



p 2 ) 



with kv > a J 2, k x > a^/2 and p G [—1,1]- The P-intensity (Af )o<t<T is specified as in ( |2.3[ ), 




i.e., we have Af = A p + Af v t + Af X t , for some A, Af , Af G R+ with A + Af + Af > 0. For 
simplicity, we assume that r t = f G R+ for all t G [0, T). Note that this specification extends the 
Heston jump-to-default model considered in [9] by allowing Af to depend on v t and on the additional 
stochastic factor X t . It can be easily checked that the specification ( |6.1[ ) satisfies As sumption [27TJ and, 
due to Prop. |2.3[ the defaultable stock price process S = (S t ) < t <T has the following dynamics: 



2.3 , the defaultable stock price process S = (S t )o<t<T nas me following d 

dS t = S t - (fi - Af ) dt + St-^Tt (p dW] + v 7 ! -P 2 dwf \ - S t - dM t 



(6.2) 



where M = (Mt) <t<r is the (P, G)-martingale defined by M t := l{r<t} — J T Af du. We also have: 



dv t = k(v — v t )dt + (jy/vt dWt 
dX t = k (x - X t )dt + o-ov^^W^ 2 



(6.3) 



6.2 Risk-neutral measures which preserve the Heston with jump-to-default 
structure 



By relying on Thm. 3.2 we now characterise the family of all risk- neutral measures Q E Q which 



preserve the Heston with jump-to-default structure, namely all risk-neutral measures Q G Q which 



leave unchanged the structure of the SDEs ( ]6.2[ )-( |673[ ) (compare also with [29], Sect. 2.4.1) 



Lemma 6.1. A risk-neutral measure Q G Q preserves the Heston with jump-to-default structure if 
and only ifdQ/dP admits the representation <\3.\^ for some ¥ -adapted processes 6 = (0 t )o< t <T and 
7 = (lt)o<t<T of the form (|3.2[) with 9 G IR 3 and 6 G IR 3x3 satisfying the following restrictions: 



(71 

r+\Q-ii-p§i 



\ 



e 



/ ©1,1 
o 



with 6i > a/2 — kv/a and 02 > <Jq/2 — k$x / o"o- 







e 



2,2 



0\ 




\ y/T- 







(6.4) 



Proof. The result follows from conditions ( |3.2| )-( |33] ) of Thm. 3.2 , noting that the preservation of the 
Heston with jump-to-default structure consists in the additional restriction Q 12 = ©2,1 = 0. □ 
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The main benefit of working with risk-neutral measures which preserve the Heston with jump-to- 
default structure consists in the possibility of obtaining closed-form solutions to the system of Riccati 
ODEs p.4[ ), as shown in the next Lemma, which follows from Lemma 10.12 of ll27l by means of 
simple (but tedious and, hence, omitted) computations. 

Lemma 6.2. Let Q G Qbe a risk-neutral measure which preserves the Heston with jump-to-default 
structure. Then the system of Riccati ODEs p.4[) admits the following solution, for all z G C 2 x zR.- 



^(t,z)=z s 

2(fcu + o-^i) 



(z 3 - zf + 2A?(1 - z 3 )) (e^* - 1) - (V2ST(e^ + l) + (a (0 M + pz 3 ) - fc) (e^* - l)) z x 

V / A7(ev / A^ t + 1) - (a (9i,i + pz 3 ) - fcj (ev^* - l) - a 2 (e^* - l) z x 
2A?(1 - z 3 )(e^* - 1) - (v^(e^ + 1) + (tr 9 2 , 2 - fc )(e^ 4 - l)) *2 



t,z) 



VA 2 (e^ + 1) - (a e 2 , 2 - k )(e^ - l) - og(eV25t _ i) Z2 



log 



/ST-(s(ei,i+p« 3 )-fc) 



/A7( eV ^Tt + i) _ (5(614 + pz 3 ) - fc)(ev / A^ t - 1) - a 2 (e^* - l)zj 



2(fc i + cro# 2 ) 



log 



/A2( e VSJt + i) _ ( CTo e 2j2 - fc )(e^ - 1) - ag(e^ - l)z 2 



+ (r + A«)(2&-l)t 



Ax := (a (0!,! + pz 3 ) - k) 2 + a 2 (z 3 - z\ + 2A?(1 - z 3 )) 
A 2 := (a o 2 , 2 - A; ) 2 + 2a 2 Af (1 - z 3 ) 



By combining the above Lemma with the results of Sect.s|4jj5} we can efficiently solve risk man- 
agement problems and compute arbitrage-free prices of general default- sensitive payoffs. 



6.3 Numerical results 



This Section reports the results of some numerical computations for the Heston with jump-to-default 
model ( |6.1[ )-( [63] ). We adopt the following parameters' specification: k = 0.565, v = 0.07, a = 0.281, 
k = 0.325, x = 0.003, a = 0.036, p = 0.1, p = —0.558. These values have been obtained in ifTOl 
by calibrating a jump-to-default model to equity options and CDS spreads on the Citigroup company 
(period: 5/2002 - 5/2006). The remaining parameters appearing in ( |6.4[ ) are specified as f = 0, 
e M = 9 2i2 = 0.002, Bi = 9 2 = 0.001 and Af = Af = A p = 0.1225. 

By relying on Proposition |4.4| and Lemma 6.2 we compute the distribution function of the default- 
able stock price Sr in the case of survival, i.e., we compute the map (T, x) H- P (Sr < x,t > T), 
for T G [0.5,3.0] and x G [0.7, 1.3], for S = I. In particular, we can observe from Fig. [Tjthat 
lim^+oo P (St < x,t > T) = P(r > T) < 1, with a more pronounced effect as T increases. 



As a second application, we construct the implied volatility surface generated by the model (6.1 ) 



( |6.3[ ). To this effect, we first compute a matrix of prices Pk(0, T) of Put options on the defaultable 
stock St, issued by a default-free third party, with maturity T G [0.5, 3.0] and moneyness K/Sq G 
[0.7, 1.3], letting A Q = 0.001 and A^ = Af , for i — 1, 2. The computation is performed via the Fast 
Fourier Transform method of Q, by relying on Corollary |5.4 and Lemma 6.2 The corresponding 
implied volatilities are then computed by using the blsimpv function in Matlab® (R2012a 64-bit 
version). Fig. [2] compares the implied volatility surface generated by the model ( |6.1[ )-( |63| ) with the 
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1, 




Figure 1: Surface (T,x) ^ P{St<x,t>T) for 
the Heston with Jump-to-Default model. 



j£ 0.28- 
^ 0.26- 




Figure 2: Implied volatility surfaces: standard Hes- 
ton (solid) and Heston + Jump-to-Default (mesh). 



implied volatility surface obtained from a standard (default-free) Heston BOB model, i.e., by letting 
\ Q = A® = A® = 0. It is evident that the introduction of a jump-to-default increases the implied 
volatility along all maturities and strikes. The increase is more pronounced for deep out-of-the-money 
options, due to the possibility of obtaining K in the case of default (compare also with eqn. ( |5.7[ )). 
There is also a strong skew effect, which tends to flatten as the maturity increases but is always more 
significant than in the default-free case. 



7 Conclusions and further developments 

We have proposed a general framework based on an affine process V and on a doubly stochastic ran- 
dom time r for the modeling of a defaultable stock. This approach allows to jointly model equity 
and credit risk, together with stochastic volatility and stochastic interest rate. Moreover, analytical 
tractability is ensured under both the physical and a risk-neutral probability measure, thanks to a char- 
acterisation of all risk-neutral measures which preserve the affine structure of (V, r). In the present 
paper, the driving process V has been specified as an affine diffusion on M.™ + x IR ' - " 1 . However, our 
approach can be easily adapted to the more general case where V is a continuous matrix- valued affine 
process (e.g., a Wishart process), as recently considered in lfT8l . 
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